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Abstract 

The influence of the boundary conditions used in the Wigner-Seitz 
(WS) approximation applied to the neutron star inner crust is exam- 
ined. The generalized energy functional method which includes the 
neutron and proton pairing correlations is used. Predictions of two 
versions of the boundary conditions are compared with each other. 
The uncertainties in the equilibrium configuration (Z,R C ) of the crust, 
where Z is the proton charge and R c , the radius of the WS cell, cor- 
respond to variation of Z by 2 - 6 units and of Rc, by 1 — 2 fm. The 
effect of the boundary conditions is enhanced at increasing density. 
These uncertainties are smaller than the variation of Z and R c com- 
ing from the pairing effects. In the case of high densities, kp > 1 fm -1 , 
the most important uncertainty occurs in the value of the neutron gap 
A n . In the WS approximation, it originates from the shell effect in 
the neutron single-particle spectrum which is rather pronounced in the 
case of larger k? and, correspondingly, small R c values, but it becomes 
negligible at lower density near the drip point. An approximate recipe 
to avoid this uncertainty is suggested. 

PACS : 26.60. +c,97.60.Jd,21.65.+f,21.60.-n,21.30.Fe 



1 Introduction 

In the last two decades the interest on the structure of the neutron star inner 
crust has been stimulated by the increasing number of observational data 
on the pulsar glitches. The latter are commonly explained in terms of the 
dynamics of superfruid vortices within the inner crust of neutron stars (see £Q 
and Refs. therein). By "inner crust" one usually indicates the part of the shell 
of a neutron star with sub-nuclear densities O.OOlpo ~ P ~ 0.5po, where po is 
the normal nuclear density. According to present-day ideas, the bulk of the 



1 



inner crust consists mainly of spherically symmetrical nuclear-like clusters 
which form a crystal matrix immersed in a sea of neutrons and virtually 
uniform sea of electrons. Such a picture was first justified microscopically in 
the classical paper by Negele and Vautherin [2] within the Wigner-Seitz (WS) 
approximation. This approximation consists in replacing the crystal by a sum 
of identical spherical cells with a nuclear cluster at the cell center. As far 
as the WS method is inconsistent with the crystal symmetry (periodicity), 
it implies a need for some artificial boundary conditions for the neutron 
wave functions at the WS cell boundary. To some extent, these conditions 
are arbitrary provided they guarantee the orthogonality and completeness of 
the single-particle basis. Unfortunately, calculations with complete taking 
into account the crystal symmetry are quite complicated. There are only 
few such calculations, being limited mainly with consideration of the deep 
(high density) layers of the crust OE], where the "lasagna" or "spaguetti" 
structure of the crust matter is assumed and the use of the WS method with 
the spherical symmetry cannot be applied. Only recently a consistent band 
theory was developed for the outer (low density) layers as well 5 J . Therefore 
up to now the WS method is quite popular and it is usually considered as the 
most practical one for systematic investigation of the inner crust structure 
in the whole density interval. For describing the matter of a neutron star 
crust, Negele and Vautherin used a version of the energy functional method 
with the density dependent effective mass m*(p). In fact, it is very close to 
the Hartree-Fock method with effective Skyrme forces. For a fixed average 
nuclear density p, the nuclear (plus electron) energy functional is minimized 
for the spherical WS cell of the radius R c . A cell contains Z protons (and 
electrons) and N=A— Z neutrons (A = (4ir/3)Rlp). In addition, the /in- 
stability condition, 

Ain - (fip + lie) = 0, (1) 

has to be fulfilled, where fi n , \i v and \i e are the chemical potentials of neutrons, 
protons and electrons, respectively. The minimization procedure is carried 
out for different values of Z and R c . The equilibrium configuration (Z, R c ) at 
the considered density corresponds to the absolute minimum in energy among 
all these possible configurations. Application of the variational principle to 
the energy functional by Negele and Vautherin for a WS cell results in a 
set of the Shrodinger-type equations for the single particle neutron functions 
(f>\(r) = R n ij(r)$ij m (n), with the standard notation. The radial functions 
Rnij( r ) obey the boundary condition at the point r = R c . As it was noted 
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above, different kinds of boundary conditions could be used. Negele and 
Vautherin used the following one: 



Rnij(r = Rc) = 



(2) 



for odd I, and 




) 



r=R c 







(3) 



for even ones. Let us denote it as BC1. The use of this version of the 
boundary conditions has been partly justified by physical considerations in 
[2], but the dependence of the results on the particular choice has never been 
discussed in detail. It is the purpose of the paper to study this problem at 
a quantitative level and to establish the corresponding uncertainty, which is 
inherent to the WS method applied to neutron star crust. For this aim, we 
compare results obtained for the BC1 with those found for an alternative kind 
of the boundary conditions (BC2) when Eq. (J2J) is valid for even I whereas 
Eq. (J3J), for odd ones. In principle, two additional kinds of the boundary 
conditions exist when Eq. J2J) or Eq. (jSJ) is used for any I. As it was noted 
in these versions have an obvious drawback for the case of the neutron 
star inner crust since they lead to an unphysical irregular behavior of the 
neutron density p n (r) in vicinity of the point r = R c . Indeed, p„(r) vanishes 
in this point in the first case and has a maximum in the second one. On 
the contrary, p n (r) is almost constant nearby the point r = R c in the case of 
the BC1 or BC2 kinds of the boundary conditions. It should be noted that 
the pairing effects were not taken into account in [2] since it was supposed 
that they are not important for the structure of the crust. The reason of 
such an assumption is the rather small contribution of the pairing effects 
to the total binding energy of the system under consideration. Recently, 
we have generalized the approach by Negele and Vautherin to describe the 
inner crust by explicitly including the neutron and proton pairing correlations 
|Hl HI El Hj in a self-consistent way. It turned out that in the whole interval of 
p the equilibrium configuration (Z, R c ) changes significantly due to pairing. 
To explain this effect, it is instructive to analyze the /3-stability condition 
Since electrons in the inner crust of a neutron star are ultra-relativistic, 
the following relation is valid: p e ~ (9nZ /4) 1 ' 3 / R c . By substituting it into 
Eq. (JTJ), one finds 



Z ~ — 



(Pn - PpfRl- 



(4) 
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The influence of pairing on the chemical potentials fi n and fi p is much 
stronger than that on the total binding energy. Their variation may be of 
the order of the gap value A ~l-2 MeV. Such a variation of fi n or \i v may 
lead to a sizable change of the equilibrium value of Z because the difference 

— is raised to the third power in Eq. (HJ). The estimate of the change 
of Z induced by this variation is as follows: 5Z=3Z5(fi n — /i p ) / — fi p ). For 
average values of k F , the difference \i n — \i v ~ 50 4- 70 MeV, hence 5Z could 
reach few units of Z. An additional change of the Z value may appear due 
to a variation of R c . Besides, as it is shown in jSJIH], the binding energy Eb 
is rather flat function of Z and different local minima Eb{Z) have often close 
values of Eb- Therefore their relative position may change after switching 
off the pairing since in general the corresponding contribution to Eb is an 
irregular function of Z. Such a situation does often occur within the WS 
approach, especially for high density values, due to the shell-type effects in 
the single-particle neutron spectrum. An example is discussed in the paper. 



2 Brief description of the method 

We used the generalized energy functional method JU] which incorporates the 
pairing effects into the original Kohn-Sham ^1] method. In this approach, 
the interaction part of the generalized energy functional depends, on equal 
footing, on the normal densities p n , p p , and the abnormal ones, z/ n , z/ p , as well: 

E int = J dr£ int (p(r), z/(r)), (5) 

where £; nt is the energy functional density. It is the sum of two components, 
the normal and the anomalous (superfluid) ones: 

£int = £norm(Pr) + £an(Ar; v t)i (6) 

where r = n,p is the isotopic index. Just as in the Kohn-Sham method, 
the prescription m* = m holds to be true. To describe the central part 
of a WS cell with the nuclear cluster inside we used the phenomenological 
nuclear energy functional £ ph by Fayans et al. JU] which describes properties 
of the terrestrial atomic nuclei with high accuracy. For describing neutron 
matter surrounding the cluster we used a microscopic energy functional £ mi 
for neutron matter based on the Argonne NN potential Vi§ [T2]. The ansatz 
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of (HI E] for the complete energy functional is a smooth matching of the 
phenomenological and the microscopic functionals at the cluster surface: 

£(p T (v),u T (r)) = £ ph (p T (r), v T {v))F m {r) + £ mi (p r (r), «r(r))(l - F n {r)\ (7) 
where the matching function F m (r) is a two-parameter Fermi function: 

F m (r) = (l + exp((r-R m )/d m ))- 1 . (8) 

Eq. (J7J) is applied both to the normal and to the anomalous components of 
the energy functional. After a detailed analysis, the matching parameters 
were chosen as follows. The diffuseness parameter was taken to be equal to 
d m =0.3 fm for any value of the average baryon density of the inner crust and 
for any configuration (Z,R C ). As to the matching radius i? m , it should be 
chosen anew in any new case, in such a way that the equality 

Pp {R m ) = 0.1p p (0) (9) 

holds. In this case, on one hand, for r < R m neutrons and protons coexist 
inside the nuclear-type cluster, and the use of a realistic phenomenological 
energy functional seems reasonable. On the other hand, at r > i? m one can 
neglect the exponentially decaying proton "tails" and consider the system 
as a pure neutron matter for which an adequate energy functional micro- 
scopically calculated can be used. The same matching parameters were used 
for normal and anomalous parts of (J7J). As far as practically all the protons 
are located inside the radius R m , the matching procedure concerns, in fact, 
only neutrons, protons being described with the pure phenomenological nu- 
clear functional. It is worth to mention that for neutron matter region, the 
ansatz is, in fact, the LDA for the microscopic part of the generalized energy 
functional. As it is commonly known, the LDA works well only provided 
the density is smoothly varying, whereas it fails in the surface region with a 
sharp density gradient. The above choice of the matching procedure and the 
values of the parameters guarantees that this region of a sharp density vari- 
ation is mainly governed by the phenomenological nuclear part of the energy 
functional which "knows how to deal with it". For the microscopic part of 
the normal component of the total energy functional (J7J) we follow to refs. 
[HI E] and take the equation of states of neutron matter calculated in 
with the Argonne Vi§ potential on the basis of Brueckner theory, taking into 
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account a small admixture of 3-body force. Its explicit form could be found 
in the cited articles. The microscopic part of the anomalous component of 
the generalized energy functional in E] was calculated for the same Vis 
potential within the Bardeen-Cooper-Schrieffer (BCS) approximation. 

3 Comparison of two kinds of boundary con- 
ditions, the BC1 versus BC2. 

In the calculations of 13 El El the boundary condition by Negele and 
Vautherin, BC1, was used. Here we repeat the analysis for the case of the 
boundary condition BC2. Results for the binding energy per a nucleon, Eb, 
are shown in Fig. 1. 

It is worth to mention that calculations for larger values of fcp should be 
considered as optional as far in this case the spaguetti phase, evidently, is 
realized. This is true for kp=1.2 fm _1 [TH E?| and, maybe, for kp=l.l fm _1 
PJ. Just as in [HlEj only even values of Z are used. The detailed comparison 
is made for A;f=0.8 fm _1 . Although the two curves Eb{Z) are quite different, 
the positions of local minima for BC1 and BC2 are close to each other, the 
distance being equal to 2 or 4 units of Z. What is of primary importance, the 
relative position of local minima for BC2 is the same as for BC1. In partic- 
ular, the positions of the absolute minimum almost coincide (Z=52 for BC1 
and Z=54 for BC2). These observations permit us to simplify calculations 
for other values of kp. In the case of BC2, we limit ourselves mainly with the 
analysis of a vicinity of the absolute minimum for BC1. The neighborhood 
of other local minima was analyzed only in the case if they have values of 
Eb{Z) close to that corresponding to the absolute minimum. It turned out 
that there is no value of kp for which the relative position of a local mini- 
mum and of the absolute one for BC1 and BC2 is different. In addition to 
systematic calculations for kp=0.6 -j- 1.2 fm" 1 , we made an extra one for a 
small density, kp=0.2 fm _1 , in vicinity of the neutron drip point. In the last 
case, two curves corresponding to BC1 and BC2 practically coincide. For all 
other values of kp the absolute minima are shifted by 2, 4 or even 6 units of 
Z. Comparison of different properties of the equilibrium configuration of the 
WS cell for various values of k-p in the case of BC1 and BC2 is presented in 
Table 1. 

There are two lines for every value of kp. The first one is given for 
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Figure 1: Binding energy per a nucleon for various k F in the BC1 case 
(solid circles connected with the solid lines) and the BC2 one (open circles 
connected with the dotted lines). 
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the Z value corresponding to the minimum of Eq in the BC1 case, the 
second one, for BC2. The only exception is /cp=0.2 fm -1 when these two 
values of Z coincide. In the last couple of columns, the average value Ap 
of the diagonal matrix element of the neutron gap at the Fermi surface is 
given. The averaging procedure involves 10 levels above /i n and 10 levels 
below. For a comparison, the infinite matter value A inf is given in the Table. 
It is calculated within the BCS approximation for the density of neutron 
matter corresponding to hp value under consideration. So drastic difference 
between A in f and Ap values in the case of /c F =0.2 fm" 1 is caused by the trivial 
reason that this kp is in the vicinity of the neutron drip point. Therefore 
the asymptotic neutron density which determines mainly the A F value is 

Table 1: Comparison of properties of equilibrium configurations of the WS 
cell for two different kinds of the boundary condition 
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significantly less than that for the uniform matter distribution. One can 
see that the influence of the boundary conditions is enhanced at increasing 
values of kp. Especially strong variation of Ap and /i n values takes place in 
the cases of kp—1.1 fm _1 and kp=1.2 fm -1 , see Fig. 2. 

To illustrate the influence of the boundary conditions on the neutron gap 
in the first case, the gap function A n (r) is drawn for both values of Z and 
both kinds of the boundary conditions. The strongest variation of the gap 
occurs in the case of Z=20. To understand the reason of such strong effect, 
we draw the neutron single particle spectrum e\ for this value of Z in Fig. 
3 for the BC1 case (the left half of the figure) and the BC2 one (the right 
one). The position of the chemical potential ji n is shown with dots. 

The two spectra are absolutely different. The reason is the shift Ae\ of 
each A-level going from BC1 to BC2. The value of this shift is approximately 




5 10 15 

r, fm 

Figure 2: The neutron gap for /cp=l-l fm -1 , Z=20 and Z=26, in the BC1 
case (solid lines) and the BC2 one (dashed lines). 
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Figure 3: The neutron single-particle spectrum e A for k F =l.l fm 1 , Z=20, 
in the BC1 case (left) and the BC2 one (right). 



equal to a half of the distance between two neighboring levels with the same 
the sign of the shift being opposite for even and odd I. The absolute 
value of the shift is proportional to and grows at increasing values of 

k F . The corresponding shifts are shown in Fig. 3 for two states, 2j 13 / 2 and 
1^23/2, which are the neighbors of ji n in the BC1 case. On average, the 
spectrum is quite dense, however in both cases there is a shell type structure 
with rather wide intervals between some neighboring levels. If one deals with 
a big inter-level space in vicinity of /i n , as in the BC1 case in Fig. 3, one 
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usually obtains a dense set of levels in this region when going to the opposite 
kind of the boundary conditions. In the BC2 case, big intervals are far from 
the Fermi surface and do not influence significantly the value of the neutron 
gap. On the contrary, in the BC1 case /i n is situated just inside such an 
interval that suppresses the gap significantly. In principle, the neutron gap 
could vanish if the interval was wider. 



1 1 1 1 I 1 1 1 1 I 1 1 1 1 I 1 1 1 1 I 1 1 1 1 I 1 




I i i i i i i i i i i ■ ■ ■ ■ i i i i i i i i i i i ■ I 
5 10 15 20 25 

r, fm 

Figure 4: The neutron gap for kp=0.8 fm -1 , Z=A2 and Z=A4, in the BC1 
case (solid lines) and the BC2 one (dashed lines). 

As to the gap function itself, for intermediate densities, kp < 1 fm -1 , 
variations caused by the choice of the boundary conditions are not significant. 
An example for the case of k-p=0.8 fm -1 is given in Fig. 4 where, just as in 
Fig. 2, four curves for A(r) are shown. The difference between any couple of 
these curves is less, of course, than the accuracy of the approach. Evidently, 
the most important uncertainty in the neutron gap value comes from using 
the BCS approximation without many-body effects like screening, which, for 
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infinite neutron matter, overestimates A n significantly jjj. 

To be definite, let us consider the "self-consistent" (Z=42 for kp=0.8 fm -1 ) 
gap function for the BC1 version of the boundary conditions as the prediction 
of the WS method for A(r) in the case of small and intermediate densities, 
kp < 1.0 fm -1 . Such a choice corresponds to that used in [9 . The "nor- 
mal" situation with the gap in the case of intermediate densities corresponds 
to much more regular single-particle neutron spectra, than that in Fig. 3. 
For the case of kp=0.8 fm -1 it is illustrated in Fig. 5. One can see that, 
although here also there are some irregularities in the inter-level distances, 
they are much smaller than the gap value A ~ 2.5 MeV. Therefore they 
don't influence the gap equation significantly. 

In the case of bigger densities, as it is seen in Fig. 2, there is a signifi- 
cant uncertainty in predictions for A(r) within the WS method. As it was 
discussed above, it originates from the shell-type structure of the neutron 
single-particle spectrum which appears in the case of high kp (small R c ) val- 
ues. We consider this effect as an artifact of the WS approximation which 
must disappear (or almost disappear) in a more advanced approach with the 
consistent consideration of the band structure. No doubts, if a forbidden 
space between two bands in vicinity of [L n in the band structure calculation 
will appear, it should be much less than the gap A value. Therefore the 
solution of the gap equation in the realistic case should be closer to that in 
the WS approximation for such a situation when there is no big inter-level 
interval close to fj, n . Returning to Fig. 2, these are the dotted line for Z=20 
and the solid one, for Z=26. Again the difference between these two curves is 
not essential, and any of them could be used as the prediction for A(r) in the 
case of fcp=l.l fm -1 . Strictly speaking, such a recipe is not a self-consistent 
one within the WS method, but it looks reasonable from the physical point 
of view. 

4 Discussion and conclusions. 

As we have seen, there are internal uncertainties inherent to the WS method 
applied to the neutron star inner crust which originate from the kind of the 
boundary condition used. In the case of very small density nearby the neu- 
tron drip point the predictions of the BC1 and BC2 versions are practically 
identical. At increasing density, with kp > 0.6 fm -1 , the uncertainty in the 
equilibrium value of Z is between 2 and 6 units, with the largest values at 
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Figure 5: The neutron single-particle spectrum E\ for /cp=0.8 fm 1 , Z=42, 
in the BC1 case (left) and the BC2 one (right). 
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the largest kp. The uncertainty in the value of R c rule, about 1 fm 

and only for kp=l.l fm -1 it turns out to be about 2 fm. However, the value 
of these uncertainties is smaller than the variations of the equilibrium config- 
uration (Z,R C ) connected with the pairing effects 9\. For the case of small 
and intermediate densities , kp < 1 fm -1 , the uncertainty in predictions 
for the gap function A(r) caused by a particular choice of the boundary 
condition is also rather small. In the case of high densities, kp > 1 fm -1 , 
the most important uncertainty occurs in the value of the neutron gap A n . 
In the WS approximation, it originates from the shell effect in the neutron 
single-particle spectrum which is rather pronounced in the case of big kp and, 
correspondingly, smaller R c values. We consider this effect as an artifact of 
the WS method which should disappear in a more consistent band structure 
calculation. We suggest an approximate recipe to avoid this uncertainty for 
the gap function A(r). We think that the most important uncertainty in 
the neutron gap value comes from the BCS approximation used. As it is 
well known, in neutron matter the many-body corrections to the BCS ap- 
proximation reduce the value of A n significantly [7j. In the next paper, we 
plan to take into account this reduction effect for more realistic description 
of neutron superfluidity in the inner crust of neutron stars. 
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for Science and Education and by the RFBR grant 06-02-17171-a. 
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